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Abstract 

Properties of disordered thin films are discussed based on the viewpoint that super- 
conducting islands are formed in the system. These lead to superconducting cor- 
relations confined in space, which are known to form spontaneously in thin films. 
Application of a perpendicular magnetic field can drive the system from the su- 
perconducting state (characterized by phase-rigidity between the sample edges) to 
an insulating state in which there are no phase-correlations between the edges of 
the system. On the insulating side the existence of superconducting islands leads to 
a non-monotonic magnetoresistance. Several other features seen in experiment are 
explained. 
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1 Introduction 



Superconductivity (SC) in disordered thin films has been a subject of intense 
study for more than a decade [Ij. Nevertheless, even elucidation of one of 
the most fundamental property of these systems, namely the superconductor- 
insulator transition (SIT), remains a puzzle. Specifically, it is still unclear 
whether this is a truely quantum phase-transition, what is the role of the mag- 
netic field in the transition, and even its universality class remains undeter- 
mined. Another profound feature which is still in debate is the non-monotonic 
magneto- resistance (MR), which in some systems can reach several orders 
of magnitude |3]. The non-monotonic MR is accompanied by a unique tem- 
perature, magnetic field and disorder dependence of the resistance PH] . On 



Preprint submitted to Elsevier 



1 February 2008 



the theoretical front, the adequacy of treating the system in terms of strictly 
bosonic excitations (so called "dirty boson" models [2]) is still questionable. 

We adopt a perspective within which the system is composed of SC islands, a 
structure implying local SC correlations. The notion of SC islands appeared 
nearly a decade ago in two contexts. Galitski and Larkin [5] suggested that 
in a strongly disordered system the SC order parameter ( and hence also Tc 
and He) fluctuate in space. As a result, at a given temperature and magnetic 
field (close to He of the clean system) there will be areas in the sample where 
the local critical field exceeds the external field. Therefore these domains still 
display SC correlations. Ghosal et.al. [B] used a locally self-consistent solution 
of the Bogoliubov-de-Gennes equations to show that even in the presence of 
extremely strong disorder there remain regions in space where the SC order 
parameter is finite, surrounded by regions of vanishingly small order parameter 
(and hence dubbed SC islands). 

As it is shown below, the concept of SC islands might explain numerous ex- 
perimental observations. In section 2 we describe the nature of the SIT as it 
emerges from the existence of SC islands with fluctuating phases. In section 
3 we show how SC islands lead to non-monotonic MR and explain its depen- 
dence on temperature and magnetic field. Section 4 is devoted to a summary 
and outlook. 



2 Prom Superconducting islands to the superconductor — insulator 
transition 

In this section we describe the nature of the SIT, based on our previous in- 
vestigations [7]. On the SC side of the transition, the system has a complex 
SC order parameter, and its amplitude strongly fluctuates in space. When a 
perpendicular magnetic field is turned on and then monotonically increased, 
it begins to penetrate the sample in the form of disordered vortices [H]. As 
a consequence, there are regions in which the SC order parameter vanishes, 
and SC islands are formed. However, these islands are still interconnected via 
a Josephson coupling (JC), and hence their phases are well correlated. This 
means that Cooper-pairs can coherently traverse through the sample and carry 
super-current, i.e. the system is in a macroscopic SC state. 

When the JC between two islands becomes smaller than the temperature, 
thermal phase fluctuations overcome the phase-locking induced by the JC, 
and the islands become separate (in terms of SC wave functions). Now, as 
the magnetic field is increased the JC between the different islands decreases. 
However, as long as there is a phase-stiff path between the edges of the sample 
(that is a path of islands with JC larger than temperature) the system will 
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maintain its SC nature. At a certain magnetic field Be, the phase-stiff path is 
broken and the system is no longer in a SC state. Since it is a highly disor- 
dered two-dimensional system, one expects that it it will become insulating, 
as indeed observed. 

An immediate consequence of the above scenario is that SC correlations still 
survive on the insulating side of the transition. This idea is corroborated by 
a number of experiments. For instance, measurements of the AC conductance 
[H] show that the superfiuid stiffness remains finite well above the transition. 
Another example is magnetic-field induced conductance oscillations with 200 
period observed in a sample with a fabricated nano-hole lattice [10] , indicating 
the presence of Cooper pairs. 

To support the above picture we have performed numerical calculations of 
phase correlations between the two edges of a highly disordered superconduct- 
ing thin film. We used a Monte-Carlo scheme [11] on-top of a self-consistent 
solution of the Bogoliubov-de-Gennes mean-field equations to obtain both the 
order parameter amplitude and phase, and calculated the correlation func- 
tion Flji = {cos{69i — 56j)) (where 56i is the shift of the order parameter 
phase from its mean-field value) as a function of magnetic field (in terms of 
the flux per plaquette, 0/0o)- Here are lattice sites on the left (right) 
edge of the sample. Details of the calculations are described elsewhere [7j. In 
Fig. [1] we show Flr (diamonds) as a function of magnetic field (or flux per 
plaquette) for a single disorder realization of a system of size 20 x 5, with 
strong disorder W/t = 1.4 (where t is the hopping integral), attractive inter- 
action strength U /t = 2, average electron density (per site) {n) = 0.92 and 
temperature T = 0.04. In addition we plot the average value of the SC order 
parameter amplitude |A| (triangles), for the same magnetic field. As can be 
seen, at a certain (disorder and density dependent) magnetic field the phase 
correlation vanish, and hence this value corresponds to B^. However, the order 
parameter amplitude is finite beyond this point and vanishes only at higher 
fields. This can also be seen from the density of states (DOS), plotted in the 
inset of Fig. [1] for two values of the magnetic field. The left inset displays the 
results for vanishing magnetic field, where a clear BCS-like DOS is observed. 
The right inset pertains to magnetic field of 0/0o = 0.06, that is above the 
SIT. However, since there are still SC correlations a pseudogap develops. This 
may indicate that our picture is relevant also to the physics of high Tc su- 
perconductors [12]. We point that similar results were obtained from tunnelig 
measurments [13], where the broadened BCS peak was accounted for by SC 
order parameter amplitude fluctuations. 

Although the transition is described in terms of thermal fluctuations, this de- 
scription does not rule out the possibility that the transition is quantum in 
nature. All one needs to do is to replace thermal fluctuations with quantum 
fluctuations and temperature with an energy scale associated with quantum 
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Fig. 1. Phase correlations between the two sides of the system Flr (diamonds) 
and the SC order parameter amphtude |A| (triangles) as a function of magnetic 
field. At a certain magnetic field Be phase-correlations vanish, indicating the SIT. 
However, |A| is still finite after this point, indicating the existence of SC islands. 
Insets : the density of states for vanishing magnetic field (left inset) and just above 
the transition (right inset), demonstrating the appearance of a "pseudogap" above 
the transition. 

fluctuations. In fact, recent measurements piT^ indicate that the transition 
crosses smoothly from a thermal to a quantum phase transition, the phe- 
nomenology of the two classes being very similar. 



3 Superconducting islands and non-monotonic magneto-resistance 

In a recent set of experiments, Sambandamurthy et.al. [3] showed that the 
magneto-resistance is non-monotonic above the SIT and develops a peak at 
a certain magnetic field -Bmax- Surprisingly, the resistance at the peak can be 
several orders of magnitude larger than the resistance at the transition (which 
is found not to be universal, see e.g. |T5|16j ). Perhaps even more surprising 
is the fact that with increasing the magnetic field beyond i?max the resistance 
drops several orders of magnitude and comes close to its value at the SIT. 
While non-monotonicity in the MR was observed more than a decade ago [17j, 
there it was a miniscule effect. Here, however, this effect is huge and cannot be 
overlooked. In addition, the resistance develops an activation dependence on 
temperature, which vanishes at low enough temperatures [3j. An experimental 
study of MR for different values of disorder [1] shows that with increasing 
disorder the value of the resistance at the MR peak increases, and at the same 
time the values of both Be and i?max decreases. 

How all these phenomena can be explained in terms of the formation of SC 
islands [18| is discussed below. The model is based on three assumptions. The 
first is that disorder induces the formation of SC islands due to fluctuations in 
the amplitude of the SC order parameter. The second assumption is that as 
the magnetic field is increased, the concentration and size of these SC islands 
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decrease. This does not mean, however, that the physical size of the islands 
(that is, the spatial extent to which the order parameter amplitude is finite) 
decreases (although it might), but rather that islands become disconnected in 
the sense of phase-fiuctuations described in the previous section. The third 
assumption is that the SC islands have a charging energy, and thus, a Cooper 
pair entering a SC island (via an Andreev tunneling process) has to overcome 
it. This charging energy is expected to be inversely proportional to the island 
size, and thus to increase with increasing magnetic field. All three assumptions 
have been corroborated by numerical calculations described elsewhere [S]. 

Consider now such a system in the strong magnetic field regime, B » 5max- 
Due to the strong magnetic field the SC islands are small and have a large 
charging energy. Now the electrons can traverse through the system in two 
types of trajectories: those which follow normal areas of the sample and do 
not cross the SC islands ("normal paths") and those in which an electron tun- 
nels into a SC island via the Andreev channel ("island paths"). The resistance 
of the normal paths, Rn-, has some value (which may depend on e.g. length, 
temperature, disorder etc.) and is assumed to be weakly affected by magnetic 
field. Due to Coulomb blockade, transport through the "island paths" is ther- 
mally activated, and hence the resistance of the island paths is of the form 
Ri ~ exp(i?c/7'), where is the typical charging energy of the island. If Ec 
is large then the main contribution to the conductance is due to transport 
along the normal paths. Consistent with experiment, the MR in this regime 
is small. 

As the magnetic field is decreased (but still in the regime B > -Bmax), niore SC 
islands are created and their size increases, but they are still small enough such 
that transport along normal paths is favourable. However, some paths which 
were normal at higher fields now become island paths and hence unavailable 
for electron transport. Thus, the effective density of normal electrons which 
contribute to the resistance diminishes, resulting in a negative MR. Eventually, 
at a certain magnetic field B = 5max some SC islands are large enough so that 
their charging energy is small and the resistance through them is comparable to 
the resistance through normal paths, i.e. Rn ~ Ri- At this point the resistance 
reaches its maximum value, since as the magnetic field is further decreased 
the SC islands are so large that transport through them is always preferred 
over transport through normal paths. An increase in number and size of the 
SC islands will thus result in a decrease in the resistance. At the critical field 
Be a phase-stiff path percolates through the system, resulting in the SIT. 

The above model was encoded into a numerical calculation using a lattice 
network model of normal resistors, SC resistors and insulating resistors (to 
mimic the Coulomb blockade) \l8l. Resistance of a normal resistor connecting 
two sites on the lattice is given by Rij = Rq exp (^^^ + 1 ^ ^ where 

Ro is a constant, rij is the distance between sites i and j, ^loc is the localization 
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length, Cj is the energy of the i-th site measured from the chemical potential 
(taken from a uniform distribution [—W/2, W/2]) and T is the temperature. 
The resistance between two (neighboring) SC sites is taken to be very small 
compared with that of a normal resistor, but still not zero and temperature 
dependent, in such a way that it vanishes as T ^ (distant SC sites are 
disconnected). The resistance between a normal site and a SC site (i.e. an 
insulating resistor) is taken to be Rns oc exp{Ec/kT), where Ec is the charging 
energy of the island. 

In Fig. [2] we plot the resistance as a function of the concentration of the SC 
islands (which corresponds to the magnetic field) for different values of temper- 
ature. The quantitative resemblance to the experimental data (inset of Fig. [2]) 
is self-evident. In addition, this simple phenomenological theory is capable of 
accounting for the breakdown of activation behaviour at low temperatures, 
the dependence of disorder and other experimental observations. 




Fig. 2. Magneto-resistance as a function of concentration of SC islands for differ- 
ent temperatures (taken from Ref. [IS]), showing good quantitative resemblance 
between the theory and experimental data (inset, taken from Ref. [3]) is evident. 



4 Summary and outlook 



In this contribution we have demonstrated that many of the properties of 
disorder SC thin films may be accounted for by the formation of SC islands. 
Specifically, we have demonstrated that phase-fluctuations between different 
islands lead to the SIT, and that above the transition formation of SC islands 
leads to a non-monotonic magneto-resistance. 

Still, there are many puzzles left in elucidating the physics of these systems. For 
instance, an unusual disorder and magnetic-field dependent anisotropy in the 
magneto-resistance was recently measured [19] and is up-to-date unexplained, 
although it seems that a percolation theory similar to that presented here may 
account for it [20]. Even more challenging is the appearance of a seemingly 
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correlated |2T] and perhaps metallic state at high magnetic field [22]. A final 
example is the fact that some materials seem to exhibit a metallic intermediate 
state at the SIT while others do not [23j. These and other questions are still 
waiting to be resolved. 



References 



[I] For a review on the superconductor-insulator transition see, e.g. A. M. Goldman 
and N. Markovic, Phys. Today 51 (11), 39 (1998). 

[2] M. P. A. Fisher, Phys. Rev. Lett. 65, 923 (1990); M. Cha, M. P. A. Fisher, S. 
M. Girvin, M. Wallin, and A. P. Young Phys. Rev. B 44, 6883 (1991). 

[3] G. Sambandamurthy, L. W. Engel, A. Johansson and D. Shahar, Phys. Rev. 
Lett. 92, 107005 (2004). 

[4] M. Steiner and A. Kapitulnik, Physica C 422, 16 (2005). 

[5] V. M. Galitski and A. I. Larkin, Phys. Rev. Lett. 87, 087001 (2001). 

[6] A. Ghosal, M. Randeria and N. Trivedi, Phys. Rev. Lett. 81, 3940 (1998); A. 
Ghosal, M. Randeria and N. Trivedi, Phys. Rev. B 65, 14501 (2001). 

[7] Y. Dubi, Y. Meir and Y. Avishai, accepted to Nature. 

[8] Y. Dubi, Ph.D. thesis (unpublished). 

[9] R. W. Crane, N. P. Armitage, A. Johansson, G. Sambandamurthy, D. Shahar, 
and G. Gruner, Phys. Rev. B 75, 094506 (2007); id, Phys. Rev. B 75, 184530 
(2007). 

[10] J. M. Valles et al., unpubhshed. 

[II] M. Mayr, G. Alvarez, C. Sen and E. Dagotto, Phys. Rev. Lett. 94, 217001 
(2005). 

[12] G. Alvarez, M. Mayr, A. Moreo and E. Dagotto, Phys. Rev. B 71, 014514 
(2005). 

[13] S.-Y. Hsu, J. A. Chernevak and J. M. Valles, Phys. Rev. Lett. 75, 132 (1995). 

[14] H. Aubin, C. A. Marrache-Kikuchi, A. Pourret, K. Behnia, L. Berge, L. 
Dumoulin and J. Lesueur, Phys. Rev. B 73, 094521 (2006). 

[15] A.F. Hebard and M.A. Paalanen, Phys. Rev. Lett. 65, 927 (1990). 

[16] N. Mason and A. Kapitulnik, Phys. Rev. B 64, 060504(R) (2001). 

[17] A.F. Hebard and M.A. Paalanen, Phys. Rev. Lett. 65, 927 (1990); D. Kowal 
and Z. Ovadyahu, Sol. St. Commun. 90, 783 (1994); V. F. Gantmakher and M. 
V. Golubkov, Pis'ma Zh. Eksp. Teor. Fiz. 61, 593 (1995) [JETP Lett. 61, 606 
(1995)]. 



7 



[18] Y. Dubi, Y. Meir and Y. Avishai, Phys. Rev. B 73, 054509 (2006). 

[19] A. Johansson, N. Stander, E. Peled, G. Sambandamurthy and D. Shahar, 
|cond-mat / 06021 60 1 (unpublished) . 

[20] Y. Meir and A. Aharony, private communications. 

[21] G. Sambandamurthy, L. W. Engel, A. Johansson, E. Peled and D. Shahar, Phys. 
Rev. Lett. 94, 017003 (2005). 

[22] Y. Yu. Butko and P. W. Adams, Nature (London) 409, 161 (2001); T. I. 
Baturina, C. Strunk, M. R. Baklanov and A. Satta, Phys. Rev. Lett. 98, 127003 
(2007). 

[23] T. I. Baturina, unpublished. 



8 



